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Abstract 

Claims concerning the possible discovery of the 0^ pentaquark, with minimal quark content 
uudds, have motivated our comprehensive study into possible pentaquark states using lattice QCD. 
We review various pentaquark interpolating fields in the literature and create a new candidate ideal 
for lattice QCD simulations. Using these interpolating fields we attempt to isolate a signal for a five- 
quark resonance. Calculations are performed using improved actions on a large 20^ x40 lattice in the 
quenched approximation. The standard lattice resonance signal of increasing attraction between 
baryon constituents for increasing quark mass is not observed for spin-^ pentaquark states. We 
conclude that evidence supporting the existence of a spin-^ pentaquark resonance does not exist 
in quenched QCD. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw 
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I. INTRODUCTION 


The recently reported observations of a baryon state with strangeness S' = +1 some 
100 MeV above the NK threshold has sparked considerable interest in excited hadron spec¬ 
troscopy. Because this state has minimal quark content uudds, its discovery would be the 
first direct evidence for baryons with an exotic quark structure — namely, baryons whose 
quantum numbers cannot be described in terms of a 3-quark conhguration alone. 


A. Phenomenology 


The experimental hndings were reported in real and quasi-real photoproduction 

experiments 0, and further positive sightings were reported in A'-nucleus collisions Ibl , vp 
and pA reactions, and in neutrino jlO, EH and deep inelastic electron scattering lH| , 
for a total of around a dozen positive results. Currently only the charge and strangeness 
of this state, which has been labeled are known; its spin, parity and isospin are as 
yet undetermined, although there are hints Hi that it may be isospin zero. The mass of 
the ©■*■ is found to be around Mq+ = 1540 MeV. However, its most striking feature is its 
exceptionally narrow width. In most cases the width has been smaller than the experimental 
resolution, while analysis of NK scattering data suggests that the width cannot be greater 
than 1 MeV Ei,E3,Ei,E3,E3. Such a narrow state, 100 MeV above threshold, presents 
a challenge to most theoretical models [l^ 000. „„„„„„ 

Subsequently, a number of null results have been reported from 
and ^ 01 colliders, as well as from pp [sH , yn j0 , hadron-p [0,0] , hadron-nucleus |34| . 
pA j35|, l36|, pA [37[, and nucleus-nucleus hxed target experiments. The production 
mechanism for the ©■*■ in these reactions would be via fragmentation, and although the frag¬ 
mentation functions are not known, these results suggest that if the ©+ exists, its production 
mechanism, along with its quantum numbers, is exotic. For more detailed accounts of the 
current experimental status of pentaquark searches see Refs. |3^ H0- 

While the experimental verihcation of the ©+ and the determination of its quantum 
numbers await dehnitive confirmation, it is timely to examine the theoretical predictions for 
the masses of S' = -|-1 pentaquark states. Numerous model studies have been carried out 
rece ntly aimed at revealing the dynamical nature of the ©"*", ranging from Sk yrm ion models 
\A‘A l43l|. QCD sum rules 
l5ll |. to name just a few. 


hadronic models 


I tne ', ranging irom aK yrin ion moaeis 
0, iH and quark models jH, il, 0,0 


B. Lattice pentaquark studies 

While models can often be helpful in obtaining a qualitative understanding of data, we 
would like to see what QCD predicts for the masses of the pentaquark states. Currently 
lattice QCD is the only quantitative method of obtaining hadronic properties directly from 
QC D.^ an^ sever^ , | imirdy ^plor atory, studies of pentaquark masses have been performed 

A crucial issue in lattice QCD analyses of excited hadrons is exactly what constitutes a 
signal for a resonance. As evidence of a resonance, most lattice studies to date have sought to 
hnd the empirical mass splitting between the ©^ and the NK threshold at the unphysically 
large quark masses used in lattice simulations. This leads to the assumption that in the 
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negative parity channel the will be about 100 MeV above the S-wave NK threshold. 
However, as we will argue below, at sufficiently large quark masses the true signal for a 
resonance on the lattice, as with all other excited states studied on the lattice jo^, E3, , 

should be the presence of binding, in which case the resonance mass would be below the 
NK threshold. 

In the positive parity channel, where the N and K must be in a relative P-wave, in 
a hnite lattice volume the energy of the NK state will typically be above the mass of 
the experimental O’'' candidate. Observation of a pentaquark mass below the P-wave NK 
threshold would then be a clear sig nal for a 0+ resonance. In all of the lattice studies, with 
the exception of Chiu & Hsieh j^, the mass of the positive parity state has been found to 
be too large to be interpreted as a candidate for the 0+. 

One should note that in obtaining a relatively low mass positive parity state, Chiu & 
Hsieh 5^ perform a linear chiral extrapolation of the pentaquark mass in using only 
the lightest few quark masses, for which the errors are relatively large. Although linear 
extrapolations of hadron properties are common in the literature, these invariably neglect 
the non-analyticities in arising from the long-range structure of hadrons associated with 


the pion cloud j65|,|66j. Csikor et al. and Sasaki use a slightly modihed extrapolation, 
in which the squared pentaquark mass is htted and extrapolated as a function of m^. As 
a cautionary note, since the chiral behaviour of pentaquark masses is at present unknown, 
extrapolation of the lattice results to physical quark masses can lead to large systematic 
uncertainties, which are generally underestimated in the literature. 

The ordering of the NK and states in the negative parity channel presents some 


challenges for lattice analyses. Sasaki and Csikor et al. argue that if the 0^ is more 
massive than the NK threshold, then one needs to extract from the correlators more than 


the lightest state with which the operators have overlap. It has been suggested that 


if one can hnd an operator that has negligible coupling to the NK state, then one can £t 
the correlation function at intermediate Euclidean times to extract the mass of the heavier 
state. The idea of simply choosing an operator that does not couple to the NK threshold 
is problematic, however, because there is no way to determine a priori the extent to which 
an operator couples to a particular state. Our approach instead will be to use a number 
of different interpolating fields, which will enhance the ability to couple to different states. 
This approach has also been adopted by Fleming |^, and by the MIT group jH^. 

Extracting multiple states in lattice QCD is usually achieved by performing a correlation 
matrix analysis, which we adopt in this work, or via Bayesian techniques. The analysis 
of Sasaki uses a single interpolating held and employs a standard analysis with an 


naiysn 

exponential £t to the correlation function. Csikor et al. , Takahashi et al. and Chiu 
& Hsieh have, on the other hand, performed correlation matrix analyses using several 
different interpolating helds. In the negative parity sector, these authors extract from their 
correlation matrices both a ground state and an excited state. In all of these studies the 
positive parity state is found to lie signihcantly higher than the negative parity ground state. 

Since the S-wave NK scattering state lies very near the lowest energy state observed on 
the lattice, the issue of extracting a genuine 0^ resonance from lattice simulations presents 
an important challenge, and a number of ideas have been proposed to distinguish between 
a true resonance state and the scattering of the free N and K states in a hnite volume 
IH^. Using the Bayesian fitting techniques, Mathur et al. js^ have examined 

the volume dependence of the residue of the ground state, noting that each state — the 
pentaquark, N, and K — is volume normalised such that the residue of the NK state is 
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proportional to the inverse spatial lattice volume. The analysis suggests that the lowest- 
lying state is the NK scattering state, but leaves open the question of the existence of a 
higher-lying pentaquark resonance state. 

However, conflicting results are reported by Alexandou et al. jH^, who hnd that ratios of 
the weights in the correlation function are more consistent with single particle states than 
scattering states. At the same time, the broad distribution of u and s quarks presented 
there suggests to us the formation of an NK scattering state. 


itte rmg 

Using a single interpolating held, Ishii et al. have introduced different boundary 

conditions in the quark propagators in an attempt to separate a genuine pentaquark reso¬ 
nance state from the NK scattering state. Here the quark propagators mix such that the 
effective mass of the NK scattering state changes, while the mass of a genuine resonance is 
unchanged. Again, the lowest-lying state displays the properties of an NK scattering state, 
but leaves open the issue of whether a higher-lying pentaquark resonance exists. 

In identifying the nature of excited states, one should also explore the possibility that 
the excited state could be a two-particle state. Since we expect that our interpolating helds 
may couple to all possible two-particle states to some degree, we compare the results of our 
correlation matrix analysis to all the possible two-particle states. In the negative parity 
sector this includes the S-wave NK, NK*, and AK* (isospin-1 only) channels, as well as 
the N'K, where N' is the lowest positive parity excitation of the nucleon. In the positive 
parity channel we consider the S-wave N*K two-particle state, where N* is the lowest-lying 
negative parity excitation of the nucleon, in addition to the P-wave NK and NK* states. 


C. Lattice resonance signatures 

Our approach to assessing the existence of a genuine pentaquark resonance is complemen¬ 
tary to the aforementioned approaches. In the following we search for evidence of attraction 
between the constituents of the pentaquark state, which is vital to the formation of a res¬ 
onance. Doing so requires careful measurement of the effective mass splitting between the 
pentaquark state and the sum of the free N and K masses measured on the same lattice. 
As discussed in detail below, attraction between the constituents of every baryon resonance 
ever calculated on the lattice Q has been sufficient to render the resonance mass 

lower than the sum of the free decay channel masses when calculated at sufficiently large 
quark mass. If the behaviour of the pentaquark is similar to that of every other resonance on 
the lattice, then searching for a signal of a pentaquark resonance above the NK threshold 
at the large quark masses typically considered in lattice QCD will mean that one is simply 
looking in the wrong place. 

One might have some concern as to whether the standard lattice resonance signature 
should appear for exotic pentaquark states where quark-antiquark annihilation cannot re¬ 
duce the quark content to a “three-quark state”. Clearly the approach to the inhnite quark 
mass limit will be different. However, the heavy quark limit is far from the quark masses 
explored in this investigation, where evidence of nontrival Fock-space components (such as 
those including qq loops) in the hadronic wave functions is abundant. For example, the 
quenched and unquenched A masses differ by more than 100 MeV at the quark masses con¬ 
sidered here with the mass lying lower in the presence of dynamical fermions. We consider 
quark masses as light as 0.05 GeV, which is much less than the hadronic scale, 1.5 GeV, 
associated with pentaquark quantum numbers. In short, the traditional resonances explored 
in lattice QGD cannot be considered simply as “three-quark states”, so that there is little 
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reason to expect the lattice resonance signature to be qualitatively different for “ordinary” 
and pentaquark baryons. 

In the process of searching for attraction it is essential to explore a large number of 
interpolating fields having the quantum numbers of the putative pentaquark state. In Sec. II 
we consider a comparatively large collection of pentaquark interpolating fields and create new 
interpolators designed to maximise the opportunity to observe attraction in the pentaquark 
state. We study two types of interpolating fields: those based on a nucleon plus kaon 
configuration, and those constructed from two diquarks coupled to an s quark. 

The technical details of the lattice simulations are discussed in Sec. Ill, where we outline 
the construction of correlation functions from interpolating fields, and the correlation matrix 
analysis, as well as the actions used in this study. It is essential to use a form of improved 
action, as large scaling violations in the standard Wilson action could lead to a false signature 
of attraction. Our simulations are therefore performed with the nonperturbatively 0(a)- 
improved FLIC fermion action ^ ^ 711, which displays nearly perfect scaling, providing 
continuum limit results at finite lattice spacing 0- The simulations are carried out on 
a large, 20^ x 40, lattic^with the (!l(a^)-tadpole-improved Luscher-Weisz plaquette plus 
rectangle gauge action [7^. The lattice spacing is 0.128 fm, as determined via the Sommer 
scale, ro = 0.49 fm, in an analysis incorporating the lattice Coulomb potential 0 

In Sec. IV we present our results for the even and odd parity pentaquark states, in both 
the isoscalar and isovector channels. As we will see, there is no evidence of attraction in the 
pentaquark channel; on the contrary, we find evidence of repulsion. As the quark masses 
increase and the quark distributions become more localised, the mass splitting between the 
lowest-lying pentaquark state and the sum of the free N and K masses is generally observed 
to increase. Moreover, the standard lattice resonance signature of the resonance mass lying 
lower than the sum of the free decay channel masses at sufficiently large quark mass is 
absent. As we conclude in Sec. V, evidence supporting the existence of a spin-i pentaquark 
resonance does not exist in quenched QCD. 


II. INTERPOLATING FIELDS 

In this section we review the interpolating fields which have been used in recent pen¬ 
taquark studies, in both the QCD sum rule approach and in lattice QCD calculations. We 
then propose new interpolators designed to maximise the opportunity to observe attraction 
between the pentaquark constituents at large quark masses. Two general types of interpo¬ 
lating fields are considered: those based on an “ViF” configuration (either nK^ or pK^), 
and those based on a “diquark-diquark-s” configuration. We examine both of these types, 
and discuss the relations between them. 


A. AA-type interpolating fields 

The simplest “AA”-type interpolating field is referred herein as the “colour-singlet” form, 

Xnk = =F {u ^ d)} , (1) 

where the =F corresponds to the isospin / = 0 and 1 channels, respectively. One can easily 
verify that the field Xnk transforms negatively under the parity transformation q —>■ Jog, 
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and therefore the negative parity state will propagate in the npper-left Dirac qnadrant of 
the correlation fnnction, contrary to the more standard “positive-parity” interpolators j^. 
Note that the colonr-index snm here corresponds to a “molecnlar” (or “fall-apart”) state 
with both the ‘W” and “ii'” components of Eq. ([Q) colonr singlets. For negative parity the 
held Xnk conples the ( large x large ) x large x ( large x large ) components of Dirac 
spinors, and shonld therefore prodnce a strong signal. For the positive parity projection (see 
Sec, nil below!. it involves one lower (or small) component, conpling 
( large x la rge ) x small x ( large x large ), which is known to lead to a weaker signal in 
this channel (62| |. 

Some anthors 


521 159| have argned that xnk is a poor choice of interpolator for accessing 


the pentaqnark resonance, and that an interpolator that snppresses the colonr-singlet NK 
channel may provide better overlap with a pentaqnark resonance, shonld it exist. Csikor et 
al. Mathnr et al. jH^, Takahashi et al. js^ and Chin et al. (in lattice calcnlations), 
and Zhn (in a QCD snm rnle calcnlation), have considered a slightly modihed form in 
which the colonr indices between the N and K components of the interpolating held are 
mixed, 
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Ta 


Cxsd ) =F {u ^ d)} 


( 2 ) 


for / = 0 and 1, respectively. We refer to this alternative colonr assignment as a “colonr- 
fnsed” interpolator, whereby the colonred three-qnark and qq pair are fnsed to form a colonr- 
singlet hadron. Of conrse, for 1/3 of the combinations the colonrs e and c will coincide, so 
that the “colonr-singlet”-“colonr-singlet” state will also arise from the held Upon 

constrncting the correlation fnnctions associated with each of these interpolators, one en- 
connters a snm of (3! x 3)^ = 324 colonr combinations with a non-trivial contribntion to 
the correlation fnnction. However, only 1/9 of these terms will be common to the colonr- 
singlet and colonr-fused correlators. It will be interesting therefore to see if increased binding 
between the pentaqnark constitnents can be observed. 

In Zhn’s QCD snm rule (QCDSR) calculation interpolating helds based on the lohe 
current were also considered, such as 


75 X 


QCDSR 


1 


= + {u ^ d) . 


It is well known that the lohe current, 

can be written as a linear combination of the standard lattice interpolator. 


( 3 ) 


( 4 ) 


( 5 ) 


and an alternate interpolator whose overlap with the ground state is suppressed by a factor 
of ~ 100 0 


^abc(^Ta^^fe)7^ _ 


( 6 ) 


In the QCD sum rule approach, the interpolator of Eq. (jHI) can be used to subtract excited 
state contributions, while the nucleon is accessed via the interpolator of Eq. ® 0 , 0 . 
However, in the lattice approach, the interpolator of Eq. Q plays little to no role in accessing 
the lowest-lying state properties for either positive or negative parities 0|. 
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B. Diquark-type interpolating fields 


The other type of interpolating field which has been discussed is one in which the non- 
strange quarks are coupled into two sets of diquarks, together with the antistrange quark. 
Jaffe and Wilczek suggested that the lowest energy diquark state would have two scalar 
diquarks in a relative P-wave coupled to the s. The lowest mass pentaquark would then be 
one containing two scalar diquarks. The configuration of two isospin 1 = 0 diquarks gives a 
purely 1 = 0 interpolating field, 

Xjw = . (7) 

If the interpolating field is local, the two diquarks in xjw are identical, but because their 
colour indices are antisymmetrised, this field vanishes identically. The field Xjw would be 
non-zero if the diquark fields were non-local. On the other hand, the use of non-local fields 
significantly increases the computational cost of lattice calculations. While this remains an 
important avenue to explore in future studies, in this work we focus on the construction of 
pentaquark operators from local fields. 

A variant of the scalar field Xjw can be utilised by observing that the antisymmetric 
tensors in Eq. o can be expanded in terms of Kronecker-5 symbols, 

^abc^dec _ ^ad^be _ ^ae^bd /g'\ 


which enables the interpolating field Xjw fo be rewritten as 

XJM/ = e“'^(n^“C'75d'){(n'^'^C'75d^)-(n^"C'75d")}C's 
One can then define two interpolating fields, 
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- 71' 


abcf^.Ta 


M'“C'75d'’)(M'"C'75d")C's 


T^Te 




( 9 ) 

( 10 ) 

( 11 ) 


which are equal but do not individually vanish. Clearly these fields transform negatively 
under parity, and, as with xnk and X^j couple ( large x large ) x ( large x large ) x large 
components for negative parity states, making them ideal for lattice simulations (62| |. 

To determine the isospin of xss, one can express the second diquark as a sum of colour 
symmetric and antisymmetric components. 


755 = C'75d'^-n^^C'75d'^) 


+ -{u^ )■ Cs 


^Te 


( 12 ) 


The first term in the braces in Eq. m, which is isoscalar, is equivalent to the field Xjw^ 
and vanishes for the reasons discussed above. The second term, which does not vanish, is 
isovector, so that the field xss is pure isospin / = 1. 

An interesting question is how much, if any, overlap exists between the diquark-type field 
Xss and the NK-tjpe fields in Sec. Iff A I This can be addressed by performing a Fierz 
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transformation on the fields. For the field Xss^ one finds: 


X55 = 




( 13 ) 


The last two terms in Eq. (USD resemble NK-tjpe interpolating helds, similar to those 
introdnced in Sec. mu while the second term corresponds to an iViF*-type conhgnration. 

Note that for the NK-\ike terms in Eq. m the colonr strnctnre corresponds to the colonr- 
singlet operator Xnk- H has been snggested that the colonr-singlet NK interpolating held 
wonld have signihcant overlap with the NK scattering state and hence not conple strongly 
to a pentaqnark resonance. On the other hand, Fierz transforming the held Xss^ analogy 
with Eq. m, wonld give rise to an NK-like term corresponding to the colonr-fnsed Xf^ 
operator. Since the helds xss and Xss eqnivalent, however, this demonstrates that 
selection of the operator x^s ("''^hh the colonr-fnsed NK conhgnration) over the operator 
Xss (with the colonr-singlet NK conhgnration) wonld be spnrions. 

Several anthors IH, |5^, have also nsed a variant of the held x.jw in lattice 
simnlations, in which a scalar diqnark is conpled to a psendoscalar diqnark jd^, 

Xps = • (14) 

In this case the two diqnarks are not identical and so the held does not vanish. Since 
both diqnarks in xps are isoscalar, this held has isospin zero and transforms positively 
nnder parity. For positive parity it conples ( large x small ) x ( large x large ) x large 
components of Dirac spinors and is therefore snitable for lattice simnlations. For the negative 
parity projection, it conples ( large x small ) x ( large x large ) x small, so that the signal 
will be donbly snppressed relative to the other negative parity state interpolators. Since it 
has been nsed in the literatnre, for completeness we also inclnde xps ia our analysis. To be 
consistent with the parity assignments of the other interpolating helds discnssed above, one 
can nse a modihed form. 


Xps 75 Xps , (15) 

which then transforms negatively nnder parity. The ehect of this is merely to switch the 
positive and negative parity components propagating in the {(1,1), (2, 2)} and {(3, 3), (4,4)} 
elements of the correlation fnnction (see Sec. mm below). 


III. LATTICE TECHNIQUES 

In this section we discnss the techniques used to extract bound state masses in lattice 
calculations. We hrst outline the computation of two-point correlation functions, both at 
the baryon level, and at the quark level in terms of the interpolating helds introduced in 
Sec. mi To study more than one interpolating held, we perform a correlation matrix analysis 
to isolate the individual states, and describe the basic steps in this analysis. Finally, we 
present some details of the lattice simulations, including the gauge and fermion actions 
used. Throughout this work we shall use the Pauli representation of the Dirac 7 -matrices 
dehned in Appendix B of Sakurai [^. In particular, the 7 -matrices are Hermitian with 
= [ 7 /., 




A. Two-point correlation fnnctions 


1. Baryon level 

Let us define a function, of the interpolating field x at Euclidean time t and momentum 


exp(-i:p • x) (0 |T x{x) x(0) | 0) , (16) 

X 

where we have suppressed the Dirac indices. Inserting a complete set of intermediate mo¬ 
mentum, energy and spin states s), 

5^exp(-ip-x) (0| x{x) \B,p',s){B,p',s\ x(0) |0) , (17) 

s,p',B X 

and using translational invariance, one can write the function Q as 

= ^exp(-EBt) ( 0 | x( 0 ) |S,p,s)(S,p,s| x( 0 ) | 0 ) , (18) 

s,B 

where Eb = \/+ fp is the energy of the state \B) and Mb is its mass. 

Despite having a specihc intrinsic parity, the interpolating held y in fact has overlap with 
both positive and negative parity states. The overlap of y with the intermediate state \B^), 
where the superscript denotes parity —1 or + 1 , at the source can be parameterised by a 
coupling strength 


(0|x(0)|5 ,p,s) = ub-{p,s) , (19) 

(0|x(0)|5+,p,s) = \j-^ 75 Ub+{p,s) , ( 20 ) 

where Ebt and Mbt correspond to the energies and masses of the negative and positive 
parity states, respectively. Note that in Eq. ( 1201 ) the matrix 75 premultiplies the spinor 
Ub+ , since the interpolating helds that we use in this analysis all transform negatively under 
parity. This is in contrast to the more standard case where the helds have positive parity 
|62| . in which case the dehnitions of in Eqs. P|) and (| 2 U|) are interchanged. 

At large Euclidean times the function Q can be written 


Q{t,^ A 


B-\B-il ■P + Mb-) , c . ^ 

A -rw-exp(—E/s-t) + A 


2 M, 




B- 


(x -p- Mb+) 
‘2Mb+ 


exp(-Es+t) , ( 21 ) 


where p is the on-shell four-momentum vector, and A^""" is the coupling strength of the held y 
at the source to the state \B^). In our lattice simulations the source is smeared, so that A'®^ 
and A^""" are not necessarily equal. With hxed boundary conditions in the time direction, 
one can project out states with dehnite parity using the matrix HE! 




1 

2 



Mb± \ 
Eb± ^7 


( 22 ) 
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The parity-projected two-point correlation function can then be expressed as the spinor 
trace of the function Q{t,p), 

G^{t,p) = tr[r^^(t,^] (23) 

= ^ exp (—• (24) 

B 

Because of the exponential suppression (at large Euclidean times) of states with energies 
greater than the ground state energy, in the large-f limit the correlation function for p = 0 
(which we use in practice in this analysis) will be dominated by the ground state with mass 
m^, 


G^(t,0) *=“ exp{-m^t) , (25) 

where Aq and Ag are the ground state couplings for the negative and positive parity states 
at the source and sink, respectively. The effective mass of the baryon is then dehned in 
terms of ratios of the correlation function at successive time slices. 




In 


G^(t -|-1,0) j 


(26) 


In the large-t limit one therefore picks out the state with the lowest mass. 


Mj(i) '=“ mi 


(27) 


In order to study masses of excited states, one can in principle attempt to £t the correlation 
function at hnite t with a sum of exponentials corresponding to the ground state plus one 
or more excited states. In practice this is difhcult, however, due to the large statistics 
required for a reliable extraction of the masses. An alternative approach is to use several 
interpolating helds, and extract masses of an orthogonal set of operators using a correlation 
matrix analysis, as we discuss in Sec. IIIIBI below. 


2. Quark level 


The two-point correlation functions discussed above are all derived at the hadronic level. 
They can be expressed in a form suitable for actual lattice simulations by substituting the 
interpolating helds in Sec. Eland performing the appropriate Wick contractions of the time- 
ordered products of helds. We use the notation (0 |Tm“(x)'U^(0)|0) = U^p{x,0) for the u 
quark, and similarly for the d and s quarks, where a and fd represent Dirac spinor indices. 

For the “molecular” NK interpolating held Xnk in Eq. m the diagonal {pK^/pK^ and 
nK^/nK^) correlation function is given by 


gNK/NK _ ^a'b'c'^d'e' gibc^de 


X < —tr 


X 


a'b'c'd'e' 


abode 


75^'^" lx, 0)75 75 /^^^ (x,0)75 


i 


tT‘='(x,0) tr 


17“'(x,0) (C'75 D'''(x,0)75C' 
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-t/- (x, 0) (^^75/^“ (x, 0)75C'j (x, 0) 

+ t/-' tr [/-'(x,0) (75 Z^'^''(x, 0)^75)^ 755 ""'*(x, 0)75 
- t/“'(x,0) (75 Z^^'’'(x, 0)75C')^ 75^""'* (:r, 0)75 (c' 75 Z^'^'(x, 0)75 )^ [/-'(x, 0) 


( 28 ) 


The propagators in Eq. (OH|) are defined from source 0 to point x, and we have used the 
relation (75S'“^(x, 0)75 )*^ = S'^“(0, x)^q, for the anti-strange quark propagator. Note that the 
first two terms in Eq. ()28j] correspond to a product of the N and K correlation functions, 
whereas the last two terms have a mixed flavour and colour structure. The correlation 
function for the operator of Eq. 0 with mixed colour labels can be obtained from 
by interchanging the c e and c' <->■ e' colour indices. 

For the pK^/nK~^ interference correlation function, one has 


-*pK^ lnK~^ _ 


^ ^ ^'b'd ^d’e' ^ 

a'b'c'd'e' 


de 


abode 


X < —tr 


17“''(x, 0)(C'75/^''“'(x, 0)^75 


ba' / 


X 17“'(x, 0)75 0)75) 0) 

+ 17^''(x,0 ) (c' 75 Z^'“'(x, 0)^75)^ t/“‘^'(x, 0)75 (75^""'*(x,0)75)^75^^"'(:r,0) 

+ 17“'(x, 0)75 (75^""'*(x, 0)75)'^ 75 D“'(x, 0)^75 (i 7 “''(x, 0))^C'757 ^'"'(x, 0) 

- 17'=''(x,0) (75/^“' (x, 0)^75)(755 ""'*(x, 0)75 ) (i 7 “^'(x, 0)75)'^C'75 /^'"'(x, 0)| , 


(29) 


for the colour-singlet interpolating field Xnk, with a similar replacement c ^ d ^ e! for 
the colour-fused field 

Similarly, the correlation function for the xss “diquark” interpolating field in Eq. (HDD is 
given by 


= - ^ ^ e“''"h'^'^C'755“'*(x,0)75C' 

a'b'c'd'e' abode 


X 1 tr 

'(f/“'(x,0))^C'75/^""'(x,0) C-f, 

tr 

'(f/“'(x,0))^C'75Zl'''(x,0) Cx5 

— tr 

'([/“'='(x, 0))'^C'75 /^“'(x, 0) ^75(17“ 

'(x,0))^C75l^""'(x,0) C75’ 


-1- tr 

'(17“^'(x,0))^C'75/^“'(x,0) ^ 75' 

tr 

'(17“'(x,0))^C'75l^"''(x,0) C-f, 

— tr 

(17““'(x,0))^C'75/^“'(x, 0) C'75(17“'(x,0))^C'75l^"''(x,0) ^75 

}■ 


(30) 

Following the parity projection in Eq. ()22|1 . the correlation functions can be made 
real by including both the U and U* gauge field configurations in the ensemble averaging 
used to construct the lattice correlation functions. This provides an improved unbiased 
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estimator which is strictly real. This is easily implemented at the correlation function level 
by observing that 


M-'{{u;}) = [cjar'au^]) (C75)“T 

holds for quark propagators. For a more detailed discussion of this issue see Refs. j63|, l7f 


B. Correlation matrix analysis 


In the previous section we described how the mass of the ground state is extracted from 
the two-point correlation function by htting a constant to the effective mass. Excited state 
masses can be extracted either by htting the correlation function with several exponentials 
(which is, in general, quite difficult to do reliably), or by using more than one interpola ting 
held. In the latter approach, which was implemented in the N* spectrum analysis in Ref. [63| 
and which we adopt in this work, a set of linearly independent operators will, in general, 
overlap with more than one state. We use a correlation matrix analysis to convert a set of 
N linearly independent operators into a set of N orthogonal operators. 

In principle, to access the entire spectrum of states would require an inhnite tower of 
operators. In practice we use a 2 x 2 correlation matrix (in particular, for the NK-tjpe 
interpolating helds), which enables us to access two states in each channel. If the analysis is 
performed at large enough Euclidean times, the contributions from the N > 2 excited states 
will be exponentially suppressed, as found in the earlier N* analysis [ej. 

Generalising the two-point correlation function in Eq. (El to the case of two diherent 
interpolating helds Xi s-iid Xj fhe sink and source, respectively, the momentum-space 
two-point correlation function matrix Gij (at p = 0) can be written as 


N-l 

= 5]] A"A"exp(-m„f) , (31) 

a=0 

where a denotes each of the N states in the tower of excited states, and we have suppressed 
the parity labels. If the operators Xj orthogonal, the matrix Gij will be diagonal, 
with the only t dependence coming from the exponential factor, in which case one would 
have a recurrence relation. 


Gijit) = exp(—rua At) Gij{t + /S.t) 5ij . (32) 

In general the operators will not be orthogonal, and a new set of operators must be cre¬ 
ated from a linear combination of the old operators using the eigenvalue equation. In the 
event that the number of states matches the number of interpolators, an orthogonal set 
of interpolators can be constructed by diagonalising the correlation matrix subject to the 
condition 


G,,(t + Af)u“ = A“G,fc(t)< , (33) 

or, 

(G-‘(«)G((+A())y = (34) 

where are real eigenvectors, and the corresponding eigenvalue is A“ = exp(—mo. At). 
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A real symmetric matrix is diagonalised by its eigenvectors. However, since our smearing 
prescriptions are different at the source and the sink, the correlation matrix is real but 
non-symmetric. Consequently, one has to solve the additional left-eigenvalue equation 

v^Gij{t + At) = X^v^Gkjit) , (35) 

for eigenvectors n", or equivalently 

< (G(i + Ai)G-‘(i))., = ( 36 ) 

The eigenvectors m" and v°‘ diagonalise the correlation matrix at times t and t + At, 

v^G,j{t + At)u^ = 

= X^vrG,j{t)u^ , (37) 

and if A" ^ X^ iov a P then, 

Gij{t + At) Uj (X 5“^ . (38) 

The projected correlation matrix v^Gij(t)uj thus describes the single state a. 

In the present analysis, for each state considered our aim will be to optimise the results at 
every quark mass. We use the covariance matrix to hnd where the y^/dof for a least squares 
£t to the effective masses is < 1.5 for all quark masses. Stepping back one time slice, we 
then apply the correlation matrix analysis. If the correlation matrix analysis is successful, 
he., the correlation matrix is invertible and the eigenvalues are real and positive, we proceed 
to the next step. If the correlation matrix analysis fails, we take another step back in time, 
and continue stepping back until the analysis is successful for a given quark mass. 

The mass of the state derived from the projected correlation matrix is then compared 
with the mass obtained using the standard analysis techniques. Any mixing of the ground 
state with excited states will result in masses from the unprojected operators which he 
between the true ground and excited state masses. Therefore, in the case of the ground 
state mass, if the new mass is smaller then we use the result derived from the correlation 
matrix; otherwise, we keep the standard analysis result. For an excited state, on the other 
hand, the result from the correlation matrix analysis is used if the new mass is larger than 
that from the standard analysis. 


C. Lattice simulations 

Having outlined the techniques used to extract excited baryon masses and the choice of 
interpolating helds, we next describe the gauge and fermion actions used in this analysis. A 
more detailed account of the actions has been given by Zanotti et al. [H^ . 


1. Gauge action 


For the gauge fields, we use the Luscher-Weisz mean-held improved plaquette plus rect¬ 


angle action [7^. The gauge action is given by 


= 


5/5 

T 


1 


^-77etr[l - f/sq(a;)] - 


P 


sq 


12^2 ^ 3 

^ rect 


77 etr[l - Urect{x )], 


(39) 
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(40a) 


where the operators Usq{x) and Urectix) are dehned as 

Usq{x) = U^{x)U,y{x + fM)Ul{x + l))Ul{x) , 

Urect{x) = U^{x) Uy{x + (j)Uy{x + v + fi) Ul{x + 2i>) ul{x + v) ul{x) 

+ f/^(a;) U^,{x + fi) Uy{x + 2/i) Ul{x + fi + v) f/^(a; + z>) Ul{x) . (40b) 

The link product Urect{x) denotes the rectangular 1x2 and 2x1 plaquettes, and for the 
tadpole improvement factor we use the plaquette measure, 

1/4 

(41) 


uo = { -7^etr(f/sq) 


The gauge conhgurations are generated using the Cabibbo-Marinari pseudoheat-bath algo¬ 
rithm with three diagonal SU(2) subgroups looped over twice. The simulations are performed 
using a parallel algorithm with appropriate link partitioning, as described in Ref. 0- 
The calculations are performed on a 20^ x 40 lattice at /3 = 4.53. The scale is set via the 
Sommer scale ro, obtained from the static quark potential FT 


V (r) = Vq + ar — e 


+ l 


where Vq, a, e and I are £t parameters, and denotes the tree-level lattice Coulomb term, 

’ll r 

- = 4vr / ^cos(k-r)Zl 44 ( 0 ,k) , (42) 

[rj J [2'kY 

with the time-time component of the gluon propagator. Note that 1144 (^ 4 , k) is 

gauge independent in the Breit frame (/C 4 = 0) since = 0. In the continuum limit, the 

Coulomb term reduces to 



(43) 


Knowledge of Vq, a and e allows one to calculate ro, which is defined by 


dV{r) 


dr 


= 1.65 


r=ro 


(44) 


Using the phenomenological value of ro = 0.49 fm, the spacing of our lattice is a = 
0.128(2) fm. 


2. Fat-link irrelevant fermion action 

For the quark helds, we use the Fat-Link Irrelevant Clover (FLICl fermion action 
which provides a new form of nonperturbative 0{a) improvement [ 8 ^. This action has 

g viously been used to study hadronic masses j^, as well as the excited baryon spectrum 
. Here fat links are generated by smearing links with their nearest transverse neighbours 
in a gauge covariant manner (APE smearing). This has the effect of reducing the problem 
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of exceptional configurations common to Wilson-style actions [^, and minimising the effect 
of renormalisation on the action improvement terms. Since only the irrelevant, higher¬ 
dimensional terms in the action are smeared, while the relevant, dimension-four operators 
are left untouched, the short-distance behaviour of the quark and gluon interactions is 
retained. The use of fat links [s^ in the irrelevant operators also removes the need to hne 
tune the clover coefficient in removing all 0(a) artifacts. 

The smearing procedure involves replacing a link, U^{x), with a sum of the link and a 
times its staples 

U,{x)^U'^{x) = {l-a)U,{x) 

4 

-|- — Ufj,{x -|- va) Ul{x + fia) + Ul{x — va) — va) U^{x — ua + pa) , 

i/=l 

V^fl 


followed by projection back to SU(3). The unitary matrix which maximises 


Tie tr 


u'; 


is selected by iterating over the three diagonal SU(2) subgroups of SU(3). The entire proce¬ 
dure of smearing followed immediately by projection is repeated n times. The fat links used 
in this work are created with a = 0.7 and n = 6, as discussed in Ref. j^. The mean-held 
improved FLIC action is given by 

, (45) 

Z.[Uq ) 

where is constructed using fat links, and the plaquette measure is calculated via 
iH) using the fat links. The factor Csw is the (Sheikholeslami-Wohlert) clover coefficient 
(84| . dehned to be 1 at tree-level. The quark hopping parameter is k = l/(2m -|- 8r), and 
we use the conventional choice of the Wilson parameter, r = 1. In Eq. iSD the mean-held 
improved Fat-Link Irrelevant Wilson action is given by 


qFL 

>JW 


Ex i’{x)ij{x) + K^ij{x) 


X,ll 


(U,{x) 


7 /. — - '<p{x + jj) 

V Uq 


UHx-fe) 


Mo 


'ipix - p) 




u, 


FL 


-'ip{x + p) + 


u 


FL 


i^ix - p) 


• (46) 


As shown by Zanotti et al. the mean-held improvement parameter for the fat links 
is very close to 1, so that the mean-held improved coefficient for Csw is adequate A 
further advantage is that one can now use highly improved dehnitions of (involving terms 

up to mJ^), which give impressive near-integer results for the topological charge |85L l8fij| . In 
particul^ we employ an O(a^) improved dehnition of as used by Bilson-Thompson et 
al. mils. 

A hxed boundary condition in the time direction is implemented by setting Ut{x,Nt) = 
0 V X in the hopping terms of the fermion action, and periodic boundary conditions are 
imposed in the spatial directions. Gauge-invariant Gaussian smearing j87| in the spatial 
dimensions is applied at the source to increase the overlap of the interpolating operators 
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with the ground states. The source-smearing technique starts with a point source, 
V’o(^o, to), at space-time location (xq, to) = (1,1,1, 8), and proceeds via the iterative scheme, 

'tpi(x,t) = '^F(x,x')'tpi_i(x',t), (47) 


1 

F{x,X ) = _j_ Q,'\ ^ ^ \j-tfii.^) ^x',x+fl + tjj^(^X /i) Sx\x—p)\^ ■ (^8) 

' ' fi=l 

Repeating the procedure N times gives the following fermion held: 

-ipNix, t) F^{x, x') iJo{x', t) . (49) 

x' 

The parameters N and a govern the size and shape of the smearing function and in our 
simulations we use N = 35 and a = 6. 

Six quark masses are used in the calculations, and the strange quark mass is taken to be 
the third largest {k = 1.2885) quark mass in each case. At this k the ss pseudoscalar mass 
is 697 MeV, which compares well with the experimental value of 2M^ — = 693 MeV 

motivated by leading order chiral perturbation theory. The smallest pion mass considered 
is = 464 MeV. We have also considered two smaller masses, but hnd that the signal 
for these becomes quite noisy, and do not include them in the analysis. The analysis for 
the VA'-type interpolators is based on a sample of 200 configurations and the analysis for 
the PS and S'S'-type interpolators is based on an ensemble of 340 conhgurations. The error 
analysis is performed by a second-order, single-elimination jackknife, with the degree 

of freedom obtained via covariance matrix hts. 


IV. RESULTS AND DISCUSSION 

In this section we present the results of our lattice simulations of pentaquark masses, in 
both the = \ and channels, and for isospin J = 0 and 1. In addition to extracting 
the masses, we also study the mass differences between the candidate pentaquark states and 
the free two-particle states. This analysis is actually crucial in determining the nature of 
the states observed on the lattice, and the identihcation of true resonances. 


A. Signatures of a resonance 

At sufficiently small quark masses, the standard lattice technology will hnd, at large 
Euclidean times, the NK decay channel as the ground state of all the hve-quark interpolating 
helds discussed in Sec. ttH In previous analyses of pentaquark masses, the observation of a 
signal at a mass slightly above the NK threshold has been interpreted j^, as a candidate 
for a pentaquark resonance. A robust test of whether a signal is a resonance or a scattering 
state should involve an analysis of the binding of the state as a function of the quark 
mass. In a simple model with an attractive potential between the meson and baryon, the 
resonance would be expected to sit lower in the potential with increasing quark mass. For 
sufficiently large quark masses, a bound state will appear {lighter than its decay products). 
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FIG. 1: Mass difference between the lowest-lying negative parity excited nucleon bound state, the 
I{J^) = ) At*(1535), and the S-wave N + tt two-particle scattering state. 



(GeV^) 

FIG. 2; Mass difference between the I{J^) = 1(1"*") A(1232) and the P-wave N + tt two-particle 
scattering state. 


and therefore become accessible using standard lattice technology. This behaviour has in 
fact been observed in every lattice study of the N* spectrum in every channel. This feature 
is central to the study of the electromagnetic properties of decuplet baryons and their 


transitions j89|, l90|, |91| in lattice QCD 


As an example, consider the = \ parity partner of the nucleon, namely the A^*(1535) 
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FIG. 3: Effective mass of the I{J^) = 0(^ ) colour singlet A^iC-type pentaquark interpolator, 
Xnk- The data correspond to k = 1.2780 (squares), 1.2885 (circles), and 1.2990 (triangles). 


baryon, in lattice QCD. In the continuum, the iV*(1535) is a resonance which decays to a 
nucleon and a pion. On the lattice, however, the A^*(| ) is stable at the (unphysically) 
larg e quark masses where its mass is smaller than the sum of the nucleon and pion masses 
[63| . To illustrate this we show in Fig. C] the mass splitting between the A^*(| ) and the 
non-interacting S-wave -|- vr two-particle state, calculated on the same lattice. For all 
values of illustrated the mass of the iV*(| ) is below that of the iV -|- tt, and the mass 
difference clearly increases in magnitude with increasing m^r. In other words, the binding 
becomes stronger at larger quark masses. 

A similar behaviour is also observed in the case of the A(1232) isobar. The mass difference 
between the = 1"^ resonance and the lowest available P-wave N + 71 two-particle energy, 
shown in Fig. |21 is negative, and, as in Fig. ^ increases with m^. In fact, this pattern is 
repeated in every other baryon channel probed in lattice QCD, such as the = \ and | 
channels, for example 6^ ^ 6^. Thus the standard lattice resonance signature for the 
resonance is the existence of a state with a mass which becomes smaller than that of the 
N + K two-particle state as the quark mass increases, with the mass difference increasing 
at larger quark masses. 


B. Negative parity isoscalar states 

We begin our discussion of the results with the isospin-0, negative parity states. The 
lowest energy of a system with a nucleon and a kaon would have the N and iP in a relative 
S-wave, in which case the overall parity is negative. Isoscalar states can be constructed with 
the iViC-type interpolating helds, xnk and X+ji^ with the PS-type held Xps- 

The effective mass for the colour singlet Xnk held is shown in Fig. |21for several repre¬ 
sentative K values. The results for the colour fused operator x+k El are very similar 
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t 

FIG. 4: As in Fig. |21 but for the I{J^) = 0(^ ) colour fused A'it'-type pentaquark interpolator, 
^NK- 



FIG. 5: As in Fig. |31 but for the I{J^) = 0(^ ) PS-type pentaquark interpolator, xps- 


to those for Xnk- (Note that to avoid clutter in the hgures we do not show the points 
at the larger t values which have larger error bars, and have little influence on the hts.) 
To determine whether these operators have signihcant overlap with more than one state, 
a correlation matrix analysis is performed. However, we hnd that it is not possible to im¬ 
prove the ground state mass as described in Sec. IIIIBl Consequently the results using the 
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TABLE I: The pion mass and the masses of the I{J^) = 0(^ ) states extracted with the colour 
singlet NK, colour fused NK and PS'-type pentaquark interpolating fields for various values of k. 


K 


aM]\fK 


aMps 

1.2780 

0 . 540 ( 2 ) 

1 . 612 ( 17 ) 

1 . 625 ( 16 ) 

1 . 601 ( 21 ) 

1.2830 

0 . 500 ( 2 ) 

1 . 539 ( 17 ) 

1 . 553 ( 16 ) 

1 . 532 ( 23 ) 

1.2885 

0 . 453 ( 2 ) 

1 . 449 ( 20 ) 

1 . 461 ( 20 ) 

1 . 458 ( 27 ) 

1.2940 

0 . 400 ( 3 ) 

1 . 349 ( 28 ) 

1 . 361 ( 29 ) 

1 . 396 ( 37 ) 

1.2990 

0 . 345 ( 3 ) 

1 . 236 ( 50 ) 

1 . 245 ( 48 ) 

1 . 372 ( 72 ) 

1.3025 

0 . 300 ( 3 ) 

1 . 145 ( 67 ) 

1 . 138 ( 80 ) 

1 . 442 ( 171 ) 


TABLE 11: Masses of the S-wave N + K, N + K*, A + K* and N' + K two-particle states. 


K 

n A/fS-wave 

VfS-wave 

Vf S-wave 

yfS-wave 

1.2780 

1 . 553 ( 10 ) 

1 . 762 ( 14 ) 

1 . 893 ( 13 ) 

2 . 180 ( 27 ) 

1.2830 

1 . 485 ( 11 ) 

1 . 704 ( 16 ) 

1 . 848 ( 14 ) 

2 . 136 ( 32 ) 

1.2885 

1 . 404 ( 13 ) 

1 . 635 ( 19 ) 

1 . 799 ( 16 ) 

2 . 092 ( 39 ) 

1.2940 

1 . 314 ( 18 ) 

1 . 561 ( 26 ) 

1 . 751 ( 19 ) 

2 . 061 ( 54 ) 

1.2990 

1 . 216 ( 32 ) 

1 . 485 ( 41 ) 

1 . 709 ( 24 ) 

2 . 056 ( 87 ) 

1.3025 

1 . 130 ( 52 ) 

1 . 421 ( 68 ) 

1 . 682 ( 29 ) 

2 . 093 ( 144 ) 


standard (Le., non-correlation matrix) analysis techniqnes are reported in this channel. For 
comparison, in Fig.^lwe also show the effective mass of the negative parity xps diqnark-type 
held. As expected, becanse of the presence of two smaller components of Dirac spinors in 
Xps compared with xnk and Xnk ’ signal here is mnch noisier than for either of the 
A^A-type helds. This is despite the fact that almost twice as many conhgnrations are nsed 
for the Xps than for the AA-type helds. 

The pentaqnark masses are extracted by htting the ehective masses in Figs. OHSl over 
appropriate t intervals, chosen according to the criterion that the x^ degree of freedom 
is less than 1.5, and preferably close to 1.0. For the Xnk and helds, htting the ehective 
mass in the window t = 19 — 23 is fonnd to optimise the y^/dof. For the xps held, the signal 
is lost at slightly earlier times, and conseqnently we ht in the time interval t = 18 — 20 . 
The results for the masses corresponding to the Xnk, Xn1< Xps helds are tabulated in 
Table m 

The pentaquark masses are compared with masses of several two-particle states, which 
are reported in Table ITTl The lowest-energy two-particle states in the I{J^) = 0(| ) channel 
are the S-wave A -|- A, the S-wave A -|- A*, the P-wave N* + A, where A* is the lowest 
negative parity excitation of the nucleon, and the S-wave A' -|- A, where A' is the hrst 
positive-parity excited state of the nucleon. The contributions to the correlation function 
from the P-wave states are likely to be suppressed, however, because there is a contribution 
to the P-wave signal from two small components of the spinors. We therefore focus on the 
S-wave states in Table El 

The positive parity A' state, which is calculated from the interpolating held 
C 1 is known to have poor overlap with the nucleon ground state, as well 
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FIG. 6: Masses of the I{J^) = 0(^ ) states extracted with the colour singlet NK^ colour fused 
NK and PS-type pentaquark interpolating fields as a function of m^. For comparison, the masses 
of the S-wave N + K, N -\- K* and N' + K two-particle states are also shown. Some of the points 
have been offset for clarity. 


TABLE III: Mass differences between the I{J^) = 0(^ ) states extracted with the NK and 
PS-type pentaquark interpolating fields and the S-wave N + K two-particle state. 


K 

uMnk - 

aMps - aMfilY" 

1.2780 

0 . 056 ( 9 ) 

0 . 051 ( 19 ) 

1.2830 

0 . 052 ( 11 ) 

0 . 052 ( 21 ) 

1.2885 

0 . 048 ( 13 ) 

0 . 060 ( 25 ) 

1.2940 

0 . 043 ( 17 ) 

0 . 086 ( 37 ) 

1.2990 

0 . 025 ( 30 ) 

0 . 148 ( 75 ) 

1.3025 

- 0 . 011 ( 67 ) 

0 . 286 ( 177 ) 


as with the low-lying I"*" excitations, snch as the Roper iV*(1440). In fact, it gives a mass 
greater than ~ 2 GeV, signihcantly above the low-lying I"*" excitations j^. We expect, 
therefore, that our pentaquark helds will not have strong overlap with the N' + K state 
either. 

The results for the extracted masses from Tableware displayed in Fig. IHlas a function 
of m^. The masses of the pentaquark states extracted from the Xnki Xnk Xps fields 
agree within the errors (the Xnk and Xfm harasses in particular are very close), although 
the errors on the xps become large at the smaller quark masses. The pentaquark masses 
are either consistent with or lie above the mass of the lowest two-particle state, namely the 
S-wave N + K. 

The mass differences AM between the low-lying pentaquark states and the two-particle 
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FIG. 7: Effective mass difference between the I{J^) = 0(^ ) state extracted with the colour 
singlet NK-iy^e pentaquark interpolator, xnk-, and the S-wave N + K two-particle state. The 
data correspond to k = 1.2780 (squares), 1.2885 (circles), and 1.2990 (triangles). 



FIG. 8: As in Fig. (3 but for the I{J^) = 0(^ ) PS'-type interpolating field, xps- 


scattering states can be better resolved by fitting the effective mass for the mass differ¬ 
ence directly. This allows for cancellation of systematic errors since the pentaquark and 
two-particle states are generated from the same gauge field configurations, and hence their 
systematic errors are strongly correlated. Figures [3 and |H1 illustrate the effective mass plots 
for the mass differences. Note that the scale of these figures is enlarged by a factor of six 
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FIG. 9: Mass difference between the I{J^) = 0(^ ) state extracted with the NK-iy^e pentaquark 
interpolating field and the S-wave N + K two-particle state. 
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FIG. 10: Mass difference between the I{J^) = 0(^ ) state extracted with the PS'-type pentaquark 
interpolating field and the S-wave N + K two-particle state. 


compared with Figs. OHSl The mass difference between the state extracted from the colour 
singlet NK interpolator and the S-wave N + K two-particle state is fitted at time slices 
t = 19 — 21, while that between the PS extracted state and the N + K state is fitted at 
t = 18 - 20. 

The results of the mass difference analysis are presented in Table irm and in Figs. M and 
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[inifor the xnk and xps fields, respectively. We see clearly that the masses derived from the 
NK pentaquark operator are consistently higher than the lowest-mass two-particle state. 
The mass difference AM is ~ 100 MeV at the larger qnark masses, and weakly dependent 
on m^, with a possible trend towards a larger AM with increasing m^. Note the size of the 
error bars for the mass difference is rednced compared with the error bars on the masses in 
Fig. El 

Since the difference between the reported experimental mass and the physical S-wave 
N + K continnnm is ~ 100 MeV, naively one may be tempted to interpret the resnlts 
in Figs. El and |H1 as a signatnre of the on the lattice. However, the behavionr of the 
pentaqnark-(V -|- K) mass difference is in marked contrast to that of all other excited N* 
states stndied on the lattice [olleiH, as discnssed in Sec. IIV Al above, for which AM is 
negative. The lack of any binding leads ns to conclnde that the observed signal is nnlikely 
to be a resonance, and may instead correspond to an iV -|- iF two-particle state. The volnme 
dependent analysis in Ref. indeed conclnded that their signal, which is consistent with 
onr results, corresponds to an NK scattering state. 


C. Negative parity isovector states 

For the isospin-1, negative parity sector, we consider three operators which can create 
I{J^) = 1(| ) states: the isovector combinations of the colour singlet Xnk and colour fused 
XXnci the S'S'-type operator, xss- As for the isoscalar case, we perform a 2 x 2 correlation 
matrix analysis for the ViF-type fields, and here we do find improved access to the lowest 
lying state. 

Using the paradigm for optimising the results described in Sec. iniBi we perform the 
correlation matrix analysis for the largest 4 quark masses starting at f = 20 with At = 3. 
Here the ground state mass is found to be lower with the correlation matrix than with the 
standard analysis, indicating that the contamination of the ground state from excited states 
is reduced. For the second lightest quark mass we fit at t = 18, and for the lightest quark 
mass at t = 17, with At = 3 in both cases. For these two lightest quark masses, the ground 
state mass is not lowered, so here the standard analysis techniques are used. For the excited 
state, the masses from the correlation matrix are all higher than with the naive analysis for 
all quark masses, thus improving the analysis. 

The effective masses for the two projected iViF-type correlation matrix states, which 
we refer to as “state 1” (for the ground state) and “state 2” (for the excited state), are 
shown in Figs. ^^and^J respectively. For comparison, we also show the effective mass plot 
for the S'S'-type field xss in Fig- ESI The ground state mass extracted with the AiF-type 
interpolator is fitted at time slices t = 22 — 26, while the mass extracted with the S'S'-type 
interpolator is fitted at time slices t = 19 — 28. 

The resulting extracted masses are tabulated in Table HVl and shown in Fig. E] A clear 
mass splitting of ~ 400 MeV is seen between the ground state and the excited state for the 
AiF-type operators. The ground state mass is consistent with that obtained from the xss 
operator for the four smallest quark masses, but is slightly smaller for the two largest quark 
masses. As for the isoscalar channel, the ground state masses are either consistent with or 
slightly above the masses of the lowest two-particle state, the S-wave A -|- A. The excited 
state lies slightly above the S-wave two-particle A -|- K* threshold, which suggests that it 
may be an admixture of A -|- K* and A -|- K* scattering states. 
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FIG. 11: Effective mass of the I{J^) = 1(5 ) state corresponding to the NK-type pentaquark 
“state 1” for several values of k = 1.2780 (squares), 1.2885 (circles), and 1.2990 (triangles). 



FIG. 12 : As in Fig. ITTl but for the I{J^) = 1 (^ ) NK-type pentaquark “state 2 ”. 


The fitted mass differences between the pentaquark and two-particle state effective masses 
are summarised in Table El where we quote the differences between the NK-type “state 
1” and the S-wave N + K, between the NK-type “state 2” and the S-wave N -|- K*, and 
between the S'S'-type and the S-wave N K. These mass differences are illustrated in 
Figs, nil nni and HU for the three cases, respectively. As for the isoscalar channel, the mass 
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FIG. 13: As in Fig. but for the I{J^) = 1(^ ) SS-type interpolator, xss- 


TABLE IV: Masses of the I{J^) = 1(| ) states extracted with the NK and 55-type pentaquark 
interpolating fields for various values of k. 


K 

aMNK(i) 

aMNK(2) 

aMss 

1.2780 

1.649(15) 

1.859(22) 

1.692(8) 

1.2830 

1.578(18) 

1.797(25) 

1.619(9) 

1.2885 

1.497(27) 

1.720(31) 

1.530(11) 

1.2940 

1.408(48) 

1.629(47) 

1.434(16) 

1.2990 

1.313(66) 

1.577(77) 

1.334(26) 

1.3025 

1.251(144) 

1.554(175) 

1.245(51) 


TABLE V: Mass differences between the I{J^) = 1(^ ) states extracted with the NK and 55- 
type pentaquark interpolating fields and the S-wave N + K, N + K* and N + K two-particle states, 
respectively. 


K 

aMNK(i) — 

aMNK{2) — 

aMss - 

1.2780 

0.067(9) 

0.109(16) 

0.100(11) 

1.2830 

0.063(13) 

0.106(18) 

0.090(15) 

1.2885 

0.059(21) 

0.099(24) 

0.071(20) 

1.2940 

0.056(35) 

0.082(39) 

0.048(26) 

1.2990 

-0.006(53) 

0.102(76) 

-0.030(78) 

1.3025 

-0.161(223) 

0.146(176) 

-0.021(75) 
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FIG. 14: Masses of the I{J^) = 1(5 ) states extracted with the NK and SS-type pentaquark 
interpolating fields as a function of compared with the masses of the S-wave N + K, N + K* 
and A + K* two-particle states. Some of the points have been offset for clarity. 


differences for the ground state are clearly positive, and weakly dependent on m^. For both 
the NK-tjpe and S'S'-type ground states, the pentaquark masses are ~ 100 MeV larger than 
the S-wave N + K two-particle state. Similarly, the difference between the excited NK-tjpe 
pentaquark and the S-wave N -|- K* is ~ 150 MeV and approximately constant with m^. 
In fact. Figs. Upland El suggest that the ~ 100 MeV mass splitting obtained at the largest 
quark mass considered is approached from below. There is no evidence of binding and no 
indication of a resonance in the I{J^) = 1(| ) channel which could be interpreted as the 
0 +. 

D. Positive parity isoscalar states 

While each of the pentaquark operators considered above transforms negatively under 
parity, they nevertheless couple to both negative and positive parity states, as discussed in 
sec.inm Here we consider whether any of the operators Xnki Xps couple to a 

bound state in the isospin-0, positive parity channel. We compare the pentaquark states 
with the masses of the lowest energy two-particle states, which correspond to the P-wave 
N + K and N -|- K*, and the S-wave N* -|- K states. 

A two-particle state in a relative P-wave can be constructed on the lattice by adding 
one unit of lattice momentum {p = 2 tt/L) to the effective mass, i?efr = \/Afgg -|- p ^, of each 
particle. This effectively raises the mass of the two-particle state relative to the positive 
parity pentaquark. If a pentaquark state exists, it should therefore clearly he below the 
lowest P-wave scattering state. 

As in the negative parity channel, we perform a correlation matrix analysis using the two 
NK-type helds in order to isolate possible excited states. While the analysis suggests the 
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FIG. 15: Mass difference between the I{J^) = 1(| ) state corresponding to the A^iiT-type pen- 
taquark “state 1” and the S-wave N + K two-particle state. 
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FIG. 16: Mass difference between the I{J^) = 1(^ ) state corresponding to the NK-iy^e pen- 
taquark “state 2” and the S-wave N + K* two-particle state. 


presence of an excited state, the signal in the positive parity channel is considerably more 
noisy than for negative parity. Consequently, in practice for this channel we revert to the 
standard analysis method and extract only the ground state. Since the colour singlet and 
colour fused operators return the same ground state mass, we present the results for the 
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FIG. 17: Mass difference between the I{J^) = 1(5 ) state extracted with the 55-type pentaquark 
interpolating field and the S-wave N + K two-particle state. 



FIG. 18; Effective mass of the I{J^) = 0(^^) colour singlet NK-type pentaquark interpolator, 
Xnk, for K = 1.2780 (squares), 1.2885 (circles) and 1.2990 (triangles). 

colour singlet operator since the signal here is less noisy. 

The effective masses for the NK and PS-type interpolators are shown in Figs. ^1 and 
fnH respectively. The signal clearly becomes noisier at earlier times, and we fit the effective 
masses for the NK-type field at f = 15 — 17, while those for the PS-type interpolator are 
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FIG. 19: As in Fig. ITHl but for the I{J^) = PS'-type pentaquark interpolator, xps- 


TABLE VI: Masses of the I{J^) = states extracted with the colour singlet NK, and 

P5-type pentaquark interpolating fields for various values of n. 


K 

aMjsiK 

aMps 

1.2780 

1.935(40) 

1.721(57) 

1.2830 

1.867(50) 

1.642(76) 

1.2885 

1.782(66) 

1.547(119) 

1.2940 

1.681(90) 

1.458(207) 

1.2990 

1.561(126) 


1.3025 

1.421(170) 



TABLE VII: The masses of the P-wave N + K, N + K* and the S-wave N* + K two-particle 


states. 


K 

Af P-wave 

VfS-wave 

VfS-wave 

1.2780 

1.692(7) 

1.891(27) 

1.873(9) 

1.2830 

1.629(8) 

1.830(32) 

1.818(9) 

1.2885 

1.558(8) 

1.760(39) 

1.755(10) 

1.2940 

1.483(10) 

1.684(53) 

1.690(11) 

1.2990 

1.414(13) 

1.594(85) 

1.631(14) 

1.3025 

1.363(17) 

1.433(134) 

1.588(17) 
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FIG. 20: Masses of the I{J^) = 0(^^) states extracted with the colour singlet NK and PS'-type 
pentaquark interpolating fields as a function of For comparison, the masses of the P-wave 
N + K and N + K* and S-wave N* + K two-particle states are also shown. Some of the points 
have been offset for clarity. 


TABLE VIII: Mass differences between the I{J^) = 0(^^) states extracted with the colour singlet 
NK and PS-type pentaquark interpolating fields and the P-wave N + K two-particle state. 


K 

aJVlNK - 

aMps - 

1.2780 

0.228(38) 

0.035(57) 

1.2830 

0.223(48) 

0.021(78) 

1.2885 

0.209(65) 

0.003(122) 

1.2940 

0.183(90) 

-0.003(210) 

1.2990 

0.132(128) 


1.3025 

0.040(174) 



fit at t = 19 — 21. 

The results are tabulated in Table EH and shown in Fig. 1201 The masses of the positive 
parity states extracted with the NK and PS-type interpolating fields are very different. 
The mass extracted with the iViF-type interpolator is similar to both the S-wave N* + K 
mass and P-wave N + K* energy, whereas the mass extracted with the PS-type interpolator 
is consistent with the P-wave N + K energy, which are given in Table IVIIIl The signal 
obtained with the PS-type interpolator is rather noisy where we fit the effective masses, 
and we therefore only present results for the four largest quark masses for this operator. As 
mentioned in Sec. IIVBI the reason the signal is so poor is that our operators do not couple 
strongly to the P-wave states due to the additional small component of the interpolating 
held spinors contributing to this state. 
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FIG. 21: Mass difference between the I{J^) = 0(^^) state extracted with the NK-iy^e pentaquark 
interpolating field and the P-wave N + K two-particle state. 
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FIG. 22: Mass difference between the I{J^) = state extracted with the PS'-type pentaquark 

interpolating field and the P-wave N + K two-particle state. 
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For the differences between the pentaquark and two-particle state masses, we also fit the 
effective masses at t = 15 — 17 for the NK-tjpe field, and t = 19 — 21 for the PS-type 
field. The results are shown in Table rvTTTl and in Figs. and |221 for the differences between 
the masses extracted with the {NK and PS-type) pentaquark interpolating fields and the 
P-wave N + K two-particle state. The mass obtained with the NK-type field is ~ 300 MeV 
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FIG. 23: Effective mass of the I{J^) = 1(^^) colour singlet NK-type pentaquark interpolator, 
Xnk- The data correspond to k = 1.2780 (squares) and 1.2885 (circles). 


heavier than the lowest energy two-particle state (P-wave N + K) for all quark masses 
considered. The mass obtained with the PS'-type held is consistent with the the lowest 
energy two-particle state (P-wave N + K) for all quark masses considered. Once again this 
suggests that there is no binding in the I{J^) = 0(i ) channel, and hence no indication of 
a 0+ resonance. 


E. Positive parity isovector states 

For the isospin-1, positive parity channel analysis, we hnd that the correlation matrix 
does not produce improved results for the ground state masses compared with the standard 
analysis. In the case of the largest three k values the algorithm requires that we step 
back three or more time slices before the correlation matrix analysis works. The use of a 
correlation matrix analysis on these data is inappropriate due to large errors in the data. 

The effective masses for the NK-tjpe and S'S'-type interpolating helds are illustrated in 
Figs. ESI and EH respectively. Because the signal for the positive parity is rather more noisy 
than in the corresponding negative parity channel, we only show the effective mass for the 
smallest and third-smallest values of k. For the NK-tjpe pentaquarks, the colour-singlet 
Xnk and colour-fused X^ fields are found to access the same ground state, and in Fig. ESI 
we only show the results of the former. 

The effective masses for the NK and S'S'-type interpolators are htted at time slices 
f = 20 — 22 for the three largest quark masses. The results for the extracted masses and the 
corresponding two-particle states are shown in Table HXl and in Fig. ESI The ground state 
masses for the NK and S'S'-type helds are again very different. The mass extracted with 
the NK-tjpe interpolator is consistent with the P-wave N + K energy, whereas the mass 
extracted with the S'S'-type interpolator is consistent with both the S-wave N* + K mass 
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FIG. 24: As in Fig. |231 but for the I{J^) = 1(^^) 5'5-type pentaquark interpolator, xss- 


TABLE IX: Masses of the I{J^) = 1(^^) states extracted with the colour singlet NK and 55-type 
pentaquark interpolating fields for various values of n. 


K 

aMjsiK 

aMss 

1.2780 

1.732(48) 

1.956(133) 

1.2830 

1.651(57) 

1.939(158) 

1.2885 

1.536(71) 

1.954(214) 


and P-wave N + K* energy. 

The results of the mass splitting analysis are shown in Table 1X1 and illustrated in Figs.^HI 
and 123 The mass difference between the NK and 55-type pentaquarks and the P-wave 
N + K two-particle state is positive for the largest quark masses. The splitting increases 
with giving an indication of repulsion, rather than attraction. In all cases, the masses 
exhibit the opposite behaviour to that which would be expected in the presence of binding. 
We therefore do not see any indication of a resonance that could be interpreted as the 
in this channel. 


F. Comparison with previous results 

To place our results in context, we summarise here the results of earlier lattice calculations 
of pentaquark masses, and compare those with the hndings of this analysis. Table lXll presents 
a concise summary of published lattice simulations, together with the results of this analysis, 
including the actions and interpolating helds used, analysis methods employed, and some 
remarks on the results. In every case, the general features of the simulation results are 
consistent with our hndings. 
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FIG. 25: Masses of the I{J^) = iCf"*") states extracted with the colour singlet NK and SS-iype 
pentaquark interpolating fields as a function of m^. For comparison, the masses of the P-wave 
N + K, S-wave N* + K and P-wave N -|- K* two-particle states are also shown. Some of the points 
have been offset for clarity. 
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TABLE X: Mass differences between the I{J^) = 1(^^) states extracted with the colour singlet 
NK and SS'-type pentaquark interpolating fields and P-wave N + K two-particle state. 


K 

aJVlNK - 

aMss - 

1.2780 

0.093(43) 

0.135(58) 

1.2830 

0.079(53) 

0.117(61) 

1.2885 

0.034(70) 

0.084(69) 


The isoscalar negative parity channel was originally presented by Csikor et al. js^ and 
Sasaki as a candidate for the G’*'. We therefore summarise in Fig. |2Hlthe results in this 
channel from the previous lattice simulations. At the larger quark masses the results of our 
analysis are in excellent agreement with those of Mathur et al. Csikor et al. and 
Ishii et al. |^. Ih^. In our analysis, and also in that of Mathur et al. improved fermion 
actions were used, and the results are in agreement at the smaller quark masses. The results 
from Csikor et al. he slightly lower than the others at small quark masses, which may 
be due to scaling violations of the Wilson fermion action. 

The central issue in all of these analyses is the interpretation of the data. The earlier work 
of Csikor et al. and Sasaki js^ identified the 0(| ) channel as a possible candidate for the 
©■*■ based on naive linear extrapolations and comparison of quenched QCD with experiment. 
Later work by Mathur et al. analysed the volume dependence of the couplings of the 
operators to this state and determined that the lowest energy state in this channel was an 
N + K scattering state. Using hybrid boundary conditions Ishii et al. also found 


35 












0.0 


0.2 0.4 0.6 

m ^ (GeV^) 


0.8 


FIG. 26: Mass difference between the I{J^) = 1(^ ) state extracted with the NK-ty^e pentaquark 
interpolating field and the P-wave N + K two-particle state. 
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FIG. 27; Mass difference between the I{J^) = state extracted with the 55-type pentaquark 

interpolating field and the P-wave N + K two-particle state. 


that this was an scattering state. Our work is consistent with the findings of both 

of these studies. 


36 












TABLE XI: Summary of published lattice QCD pentaquark studies, including the helds used, a brief description of the analysis techniques, 
and some observations from the work. 


Group 

Action 

Operators 

Analysis methods 

Observations 

Lasscock et al. 

FLIC 

X.NK 1 

XSSy XPS 

2x2 correlation matrix; 
mass splittings analysis with 
NK, NK* 

I{J^) = 1 ( 5 *) NK scattering states; 

I{J^) = 0(^ ) NK scattering state; 

0(|^) N*K scattering state 

Csikor et al. 521 

Wilson 

Xnk 1 Xx^ 

2x2 correlation matrix; 

mass ratio with NK 

0(^ ) NK degenerate state; 

0(^ ) excited state, deemed too massive 

Sasaki 

Wilson 

XPS 

standard analysis 

0(^ ) above S-wave NK] 

0(^^) above P-wave NK 

Mathur et al. 541 

overlap 

Xnk 1 Xfj^ 

volume dependence 

0,1(^ ) NK scattering state; 

0,P-wave NK degenerate 

Ishii et al. 

S3, 50] 

Wilson 

XPS 

hybrid boundary conditions; 
Bayesian analysis 

0(^ ) NK scattering state; 

0(|^) deemed too massive 

Alexandrou et al. 

Wilson 

XPS 

volume dependence 

0(| ) more consistent with single particle state; 
NK scattering state not seen 

Chiu et al. 

58] 

domain wall 

Xnk •) 

XPS 

3x3 correlation matrix 

0(^ ) NK scattering state; 
ground state 0(^ ) less massive than 

Takahashi et al. 'SS] 

Wilson 

Xnk 1 Xxix 

2x2 correlation matrix; 
volume dependence 

i NK scattering state; 

^ excited state; 

N*K scattering state 


“The AA-type fields used by Chiu et al. differ by a 75 in the nucleon part of the operator from the other groups listed, which effectively reverses 
the intrinsic parity of the operator. 
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FIG. 28; Compilation of results for the lowest-lying I{J^) = 0(^ ) state from lattice QCD 
pentaquark studies. 


V. CONCLUSION 

We have performed a comprehensive analysis of interpolating helds holding the promise 
to provide good overlap between the the QCD vacnnm and low-lying pentaqnark states. 
Central to onr analysis is the search for evidence of attraction between the constitnents of 
pentaqnark states as the inpnt qnark masses increase. Every other baryon resonance ever 
stndied on the lattice becomes stable on the lattice at snfficiently large qnark masses. This 
is the standard resonance signatnre in lattice QCD. The mass of the resonance becomes less 
than the snm of the masses of its decay prodncts and is prevented from decaying by energy 
conservation. Attraction is essential to the formation of a resonance in the light qnark mass 
regime of QCD. 

Onr resnlts reveal no evidence of attraction that leads to a bonnd pentaqnark state at large 
qnark masses. Rather, evidence of repnlsion is evident in the correlation fnnctions giving 
rise to the lowest-lying pentaqnark masses. This is particnlarly evident in the I{J^) = 1(| ) 
state and in the more accnrate resnlts for the 0(i ) state. Similarly, both positive parity 
states show an increasing mass splitting between pentaqnark and two-particle states, again 
snggesting repnlsion as opposed to attraction. 

Moreover, in every case where an interpolating held was constructed to favor 
states which are more exotic than the colour-singlet paring of a iC and N, the approach to 
the lowest-lying state was compromised. In most cases, the same ground state mass was 
recovered in the correlation function analysis, but with increased error bars. This provides 
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further evidence that the lowest lying state is simply an scattering state. 

In the case of the I{J^) = 1(| ) state, the colour-fused NK interpolator of Eq. 0 
had sufficient overlap with an exited state to allow a successful correlation matrix analysis. 
Again, the exotic colour-fused NK interpolating failed to produce evidence of a bound 
pentaquark state, the signature of a resonance on the lattice. 

Similarly, the scalar diquark-type interpolating held of Eq. m produced effective masses 
that he higher than those recovered from the colour-singlet NK-tjpe interpolating held of 
Eq. (HI). Again, a low-lying pentaquark state was not accessed, indicating the absence of the 
standard lattice resonance signature. In short, evidence supporting the existence of a spin-| 
pentaquark resonance does not exist in quenched QCD. 

This result makes it clear that a similar analysis in full dynamical-fermion QCD is es¬ 
sential to resolving the fate of the putative pentaquark resonance. We have resolved mass 
splittings of the order of 100 MeV, and one might wonder what ehect the dynamics of 
full QCD could have on this state. As diherences in self-energies between full and quenched 
QCD of order 100 MeV or more have been observed one cannot yet rule out the possible 
existence of a pentaquark in full QCD. 
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